These summaries are based on data at the block group level from the years 1980-2020. Data are also available at the block, tract, and county levels.
meta %>%
filter(su_tract == 1) %>%
select(varname, about) %>% as.list()
## $varname
## [1] "locationname" "countyname"
## [3] "totalpopulation" "Coronary_Heart_Disease2018"
## [5] "Binge_Drinking2018" "Mental_Health2018"
## [7] "High_Blood_Pressure2017" "Physical_Inactivity2018"
## [9] "Diabetes2018" "Current_Smoking2018"
## [11] "Cancer_except_skin2018" "Current_Asthma2018"
## [13] "Dental_Visit2018" "High_Cholesterol2017"
## [15] "COPD2018" "Obesity2018"
## [17] "Physical_Health2018" "Health_Insurance2018"
## [19] "less_than_sevenhr_sleep2018" "Annual_Checkup2018"
##
## $about
## [1] "Tract FIPS code"
## [2] "County name"
## [3] "Total population of the tract as of the 2010 census"
## [4] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had angina or coronary heart disease. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [5] "Adjusted percent of respondents aged >= 18 years who report having five or more drinks (men) or four or more drinks (women) on an occassion in the past 30 days. Not an indicator of the frequency of binge drinking or the specific amount of alcohol consumed."
## [6] "Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their mental health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
## [7] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have high blood pressure. Women who were told high blood pressure only during pregnancy and those who wree told they had borderline hypertension were not included. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate. Not a measure of current high blood pressure."
## [8] "Adjusted percent of respondents aged >= 18 years who answered \"no\" to the following question: \"During the past month, other than your regular job, did you participate in any physical activity or exercise such as running, calisthenics, golf, gardening, or walking for exercise?\" Only captures information about nonoccupational physical activity."
## [9] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had diabetes (other than diabetes during pregnancy). Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [10] "Adjusted percent of respondents aged >= 18 years who report having smoked >= 100 cigarettes in their lifetime and currently smoke every day or some days. Not a measure of lifetime or current number of cigarettes smoked, and each of these factors can affect the risk for acquiring chronic disease from smoking. Not not measure intent or attempts to quit smoking among smokers or exposure to secondhand smoke among nonsmokers."
## [11] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have any type of cancer except skin cancer. Not specific to cancer type. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [12] "Adjusted percent of respondents aged >= 18 years who answer \"yes\" to both of the following questions: (1) \"Have you ever been told by a doctor, nurse, or other health professional that you have asthma?\" and (2) \"Do you still have asthma?\" This indicator requires doctor diagnosis, which may not include all persons with asthma"
## [13] "Adjusted percent of respondents aged >= 18 years who report having been to the dentist or dental clinic in the previous year."
## [14] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had high cholesterol out of all respondents who have had their cholesterol checked in the last 5 years. Reliability can be low on this indicator because patients might not know what specific tests have been performed on their blood samples unless their is an issue and many patients cannot afford to have their cholesterol checked."
## [15] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had COPD, emphysema, or chronic bronchitis. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [16] "Adjusted percent of respondents aged >= 18 years who have a BMI >= 30 kg/m^2 calculated from self-reported weight and height excluding respondents who were <3ft tall or >= 8ft; weighed <50lbs or >= 650 lbs; BMI < 12 or >= 100; pregnant women. Self-reports of height and weight lead to lower BMI estimates compared to height and weight measurements."
## [17] "Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their physical health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
## [18] "Adjusted percent of respondents aged 18-64 who report having no current health insurance coverage. All persons >= 65 are eligible for Medicare. This indicator is likely to be an underestimate due to how variable health insurance coverage can be."
## [19] "Adjusted percent of respondents aged >= 18 years who report usually getting insufficient sleep (<7 hours for those aged >= 18 years, on average, during a 24-hour period). Indicator does not measure variations in sleep duration, quality of sleep, or specific sleep problems."
## [20] "Adjusted percent of respondents aged >= 18 years who report having been to a doctor for a routine checkup (e.g., a general physical exam, not an exam for specifiic injury, illness, condition) in the previous year."
glimpse(cvldat)
## Rows: 50
## Columns: 20
## $ locationname <dbl> 51540000600, 51540000800, 51540001000, 515…
## $ countyname <chr> "Charlottesville", "Charlottesville", "Cha…
## $ totalpopulation <int> 4351, 3642, 2783, 3324, 3305, 2791, 4278, …
## $ Coronary_Heart_Disease2018 <dbl> 2.2, 7.1, 5.1, 5.8, 5.9, 3.5, 3.9, 3.9, 3.…
## $ Binge_Drinking2018 <dbl> 19.2, 15.8, 19.1, 16.7, 16.2, 21.2, 20.7, …
## $ Mental_Health2018 <dbl> 21.1, 11.2, 10.3, 16.2, 18.2, 15.0, 12.9, …
## $ High_Blood_Pressure2017 <dbl> 13.6, 33.7, 27.8, 33.5, 33.4, 21.6, 20.6, …
## $ Physical_Inactivity2018 <dbl> 23.1, 21.2, 16.1, 27.2, 29.8, 18.8, 16.2, …
## $ Diabetes2018 <dbl> 4.0, 10.5, 7.9, 12.1, 12.8, 6.7, 5.9, 6.6,…
## $ Current_Smoking2018 <dbl> 19.1, 14.2, 11.5, 21.8, 24.5, 15.4, 12.2, …
## $ Cancer_except_skin2018 <dbl> 1.4, 8.4, 6.9, 5.0, 4.4, 3.6, 4.7, 5.4, 1.…
## $ Current_Asthma2018 <dbl> 10.4, 8.7, 8.0, 10.7, 11.3, 9.5, 8.8, 8.1,…
## $ Dental_Visit2018 <dbl> 56.6, 73.2, 77.6, 60.0, 54.1, 68.9, 75.6, …
## $ High_Cholesterol2017 <dbl> 14.9, 36.1, 31.4, 31.0, 30.1, 23.3, 25.6, …
## $ COPD2018 <dbl> 4.0, 6.7, 4.5, 7.0, 7.6, 4.2, 4.1, 4.0, 5.…
## $ Obesity2018 <dbl> 24.1, 27.6, 26.2, 35.9, 37.7, 26.5, 23.2, …
## $ Physical_Health2018 <dbl> 9.6, 11.3, 8.7, 13.8, 15.2, 8.9, 8.0, 7.8,…
## $ Health_Insurance2018 <dbl> 19.4, 12.7, 9.8, 18.6, 21.1, 14.1, 10.8, 1…
## $ less_than_sevenhr_sleep2018 <dbl> 40.2, 34.4, 33.1, 42.3, 44.3, 36.9, 33.7, …
## $ Annual_Checkup2018 <dbl> 72.1, 80.5, 78.1, 78.2, 77.7, 74.7, 75.8, …
cvldat %>% select(Coronary_Heart_Disease2018, Binge_Drinking2018, Mental_Health2018, High_Blood_Pressure2017, Physical_Inactivity2018, Diabetes2018, Current_Smoking2018,
Cancer_except_skin2018, Current_Asthma2018, Dental_Visit2018, High_Cholesterol2017, COPD2018, Obesity2018, Physical_Health2018, Health_Insurance2018, less_than_sevenhr_sleep2018,
Annual_Checkup2018) %>%
select(where(~is.numeric(.x))) %>%
as.data.frame() %>%
stargazer(., type = "text", title = "Summary Statistics", digits = 2,
summary.stat = c("mean", "sd", "min", "median", "max"))
##
## Summary Statistics
## =============================================================
## Statistic Mean St. Dev. Min Median Max
## -------------------------------------------------------------
## Coronary_Heart_Disease2018 5.92 1.78 1.30 6.00 9.00
## Binge_Drinking2018 17.22 1.89 14.40 16.75 21.50
## Mental_Health2018 13.55 3.92 8.30 12.80 29.70
## High_Blood_Pressure2017 30.03 6.89 11.00 31.75 38.50
## Physical_Inactivity2018 21.54 4.59 14.00 21.25 31.60
## Diabetes2018 9.36 2.73 3 9.4 14
## Current_Smoking2018 16.52 4.37 8.70 15.25 24.50
## Cancer_except_skin2018 6.58 2.11 0.80 7.20 9.20
## Current_Asthma2018 9.28 1.14 7.40 9.15 13.20
## Dental_Visit2018 70.31 7.92 46.20 71.30 83.70
## High_Cholesterol2017 32.01 6.92 10.70 34.25 39.80
## COPD2018 6.43 1.80 3 6.2 11
## Obesity2018 30.14 5.16 20.10 28.40 40.60
## Physical_Health2018 11.49 2.42 7.10 11.35 16.80
## Health_Insurance2018 13.72 3.78 6.70 13.55 23.80
## less_than_sevenhr_sleep2018 35.75 3.45 28.70 35.65 44.30
## Annual_Checkup2018 78.33 2.67 71.50 79.05 82.10
## -------------------------------------------------------------
These distributions show data at the block group level from the years 1980-2020. Data are also available at the block, tract, and county levels.
longdat <- cvldat %>% select(c(locationname, Coronary_Heart_Disease2018, Binge_Drinking2018, Mental_Health2018, High_Blood_Pressure2017, Physical_Inactivity2018, Diabetes2018, Current_Smoking2018,
Cancer_except_skin2018, Current_Asthma2018, Dental_Visit2018, High_Cholesterol2017, COPD2018, Obesity2018, Physical_Health2018, Health_Insurance2018, less_than_sevenhr_sleep2018,
Annual_Checkup2018)) %>% pivot_longer(-locationname, names_to = "measure", values_to = "value")
longdat$measure <- factor(longdat$measure,
levels = c("Coronary_Heart_Disease2018", "Binge_Drinking2018", "Mental_Health2018", "High_Blood_Pressure2017", "Physical_Inactivity2018", "Diabetes2018", "Current_Smoking2018", "Cancer_except_skin2018", "Current_Asthma2018", "Dental_Visit2018", "High_Cholesterol2017", "COPD2018", "Obesity2018", "Physical_Health2018", "Health_Insurance2018", "less_than_sevenhr_sleep2018", "Annual_Checkup2018"))
longdat %>%
ggplot(aes(x = value, fill = measure)) +
scale_fill_viridis(option = "Blues", discrete = TRUE, guide = FALSE) +
geom_histogram() +
facet_wrap(~measure, scales = "free", ncol = 4)
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning in viridisLite::viridis(n, alpha, begin, end, direction, option): Option
## 'Blues' does not exist. Defaulting to 'viridis'.
meta %>%
filter(varname %in% c("Coronary_Heart_Disease2018", "Binge_Drinking2018", "Mental_Health2018", "High_Blood_Pressure2017", "Physical_Inactivity2018", "Diabetes2018", "Current_Smoking2018", "Cancer_except_skin2018", "Current_Asthma2018", "Dental_Visit2018", "High_Cholesterol2017", "COPD2018", "Obesity2018", "Physical_Health2018", "Health_Insurance2018", "less_than_sevenhr_sleep2018", "Annual_Checkup2018")) %>%
mutate(label = paste0(varname, ": ", about)) %>%
select(label) %>%
as.list()
$label [1] "Coronary_Heart_Disease2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had angina or coronary heart disease. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[2] "Binge_Drinking2018: Adjusted percent of respondents aged >= 18 years who report having five or more drinks (men) or four or more drinks (women) on an occassion in the past 30 days. Not an indicator of the frequency of binge drinking or the specific amount of alcohol consumed."
[3] "Mental_Health2018: Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their mental health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
[4] "High_Blood_Pressure2017: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have high blood pressure. Women who were told high blood pressure only during pregnancy and those who wree told they had borderline hypertension were not included. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate. Not a measure of current high blood pressure." [5] "Physical_Inactivity2018: Adjusted percent of respondents aged >= 18 years who answered "no" to the following question: "During the past month, other than your regular job, did you participate in any physical activity or exercise such as running, calisthenics, golf, gardening, or walking for exercise?" Only captures information about nonoccupational physical activity."
[6] "Diabetes2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had diabetes (other than diabetes during pregnancy). Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[7] "Current_Smoking2018: Adjusted percent of respondents aged >= 18 years who report having smoked >= 100 cigarettes in their lifetime and currently smoke every day or some days. Not a measure of lifetime or current number of cigarettes smoked, and each of these factors can affect the risk for acquiring chronic disease from smoking. Not not measure intent or attempts to quit smoking among smokers or exposure to secondhand smoke among nonsmokers."
[8] "Cancer_except_skin2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have any type of cancer except skin cancer. Not specific to cancer type. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[9] "Current_Asthma2018: Adjusted percent of respondents aged >= 18 years who answer "yes" to both of the following questions: (1) "Have you ever been told by a doctor, nurse, or other health professional that you have asthma?" and (2) "Do you still have asthma?" This indicator requires doctor diagnosis, which may not include all persons with asthma"
[10] "Dental_Visit2018: Adjusted percent of respondents aged >= 18 years who report having been to the dentist or dental clinic in the previous year."
[11] "High_Cholesterol2017: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had high cholesterol out of all respondents who have had their cholesterol checked in the last 5 years. Reliability can be low on this indicator because patients might not know what specific tests have been performed on their blood samples unless their is an issue and many patients cannot afford to have their cholesterol checked."
[12] "COPD2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had COPD, emphysema, or chronic bronchitis. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[13] "Obesity2018: Adjusted percent of respondents aged >= 18 years who have a BMI >= 30 kg/m^2 calculated from self-reported weight and height excluding respondents who were <3ft tall or >= 8ft; weighed <50lbs or >= 650 lbs; BMI < 12 or >= 100; pregnant women. Self-reports of height and weight lead to lower BMI estimates compared to height and weight measurements."
[14] "Physical_Health2018: Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their physical health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
[15] "Health_Insurance2018: Adjusted percent of respondents aged 18-64 who report having no current health insurance coverage. All persons >= 65 are eligible for Medicare. This indicator is likely to be an underestimate due to how variable health insurance coverage can be."
[16] "less_than_sevenhr_sleep2018: Adjusted percent of respondents aged >= 18 years who report usually getting insufficient sleep (<7 hours for those aged >= 18 years, on average, during a 24-hour period). Indicator does not measure variations in sleep duration, quality of sleep, or specific sleep problems."
[17] "Annual_Checkup2018: Adjusted percent of respondents aged >= 18 years who report having been to a doctor for a routine checkup (e.g., a general physical exam, not an exam for specifiic injury, illness, condition) in the previous year."
pal <- colorNumeric("Blues", domain = mapdat$Current_Asthma2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Current_Asthma2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with asthma: ", mapdat$Current_Asthma2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Current_Asthma2018),
title = "Predicted % of adults <br> with current asthma <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Coronary_Heart_Disease2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Coronary_Heart_Disease2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with Coronary Heart Disease: ", mapdat$Coronary_Heart_Disease2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Coronary_Heart_Disease2018),
title = "Predicted % of adults <br> with Coronary Heart Disease <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Cancer_except_skin2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Cancer_except_skin2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with cancer: ", mapdat$Coronary_Heart_Disease2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Cancer_except_skin2018),
title = "Predicted % of adults <br> with cancer (excluding skin) <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Diabetes2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Diabetes2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with diabetes: ", mapdat$Diabetes2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Diabetes2018),
title = "Predicted % of adults <br> with diabetes <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Health_Insurance2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Health_Insurance2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults without health insurance: ", mapdat$Health_Insurance2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Health_Insurance2018),
title = "Predicted % of adults <br> without health insurance <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Annual_Checkup2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Annual_Checkup2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who went <br> to the doctor in last year: ", mapdat$Annual_Checkup2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Annual_Checkup2018),
title = "Predicted % of adults <br> who went to the doctor <br> in the last year <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Dental_Visit2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Dental_Visit2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who went <br> to the dentist in last year: ", mapdat$Dental_Visit2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Dental_Visit2018),
title = "Predicted % of adults <br> who went to the dentist <br> in the last year <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Current_Smoking2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Current_Smoking2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who smoke: ", mapdat$Current_Smoking2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Current_Smoking2018),
title = "Predicted % of adults <br> who currently smoke <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$Obesity2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Obesity2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who have a BMI >= 30: ", mapdat$Obesity2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Obesity2018),
title = "Predicted % of adults <br> who had a BMI >= 30 <br> in 2018", opacity = 0.7)
pal <- colorNumeric("Blues", domain = mapdat$less_than_sevenhr_sleep2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$less_than_sevenhr_sleep2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who average <7 hrs <br> of sleep a night in 2018: ", mapdat$less_than_sevenhr_sleep2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$less_than_sevenhr_sleep2018),
title = "Predicted % of adults <br> who average <7 hrs <br> of sleep a night in 2018", opacity = 0.7)